Finite size effects and the order of a phase transition in fragmenting nuclear systems 
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We discuss the implications ol finite size effects on the 
determination of the order of a phase transition which may 
occur in infinite systems. We introduce a specific model to 
which we apply different tests. They are aimed to characterise 
the smoothed transition observed in a finite system. We show 
that the microcanonical ensemble may be a useful framework 
for the determination of the nature of such transitions. 
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Since the first attempt to use fragmentation experi- 
ments on nuclei pl-Q in order to construct the caloric 
curve which links the temperature to the excitation en- 
ergy, many efforts have been made by different groups in 
order to put the experimental results on firm grounds. 
Further work is certainly necessary in order to eliminate 
as much as possible the remaining uncertainties which 
are present and difficult to get under control. At the 
present stage it is not possible to see a clear-cut sign for 
the existence of a phase transition and hence even more 
difficult to get information about its order. 

The approximate plateau observed in the first exper- 
imental measurements was interpreted as a sign for the 
presence of a first order liquid-gas phase transition [y. 
This was in agreement with the theoretical expectation 
that the equation of state should present the character- 
istics of a liquid-gas transition j^. However, the latest 
analyses obtained by means of peripheral collisions show 
a monotonously increasing curve and the onset of a steep 
rise on the high energy side of the caloric curve |^ . 

This last result [|[0"|i] ^a-y be interpreted as a sign for 
a second order phase transition. It cannot be explained 
with the most commonly used theoretical models showing 
a liquid-gas transition (Fisher droplet model ||l^, lattice 
gas model (LGM) uM), which present a single critical 
point in the phase diagram. This raises the question of 
the adequacy of the LGM which is the simplest model 
describing a short-ranged interaction between classical 
particles inside a fixed volume. A possible explanation 
of this fact has been proposed recently by the authors 



of Ref . [E2| . They observed that finite size effects could 
produce a scaling behaviour in fragment observables in- 
side the coexistence region. These signals appear to be 
suppressed in the thermodynamic limit. If this scaling 
behaviour is related to the thermodynamic transition it 
shows that a finite system can present misleading indica- 
tions concerning the order of the transition in the infinite 
system. However, these indications are physical when one 
deals with small systems like fragmenting nuclei. In this 
sense one can speak about "phase transitions in finite 
systems" or "crossovers" (smoothed transitions). In the 
following we shall elaborate on this subject. 

The LGM is usually formulated in the framework of 
the grand canonical ensemble. The basic variables of 
this model are the temperature T and the density p of 
particles. In the grand canonical formulation, the num- 
ber of particles is not strictly fixed, but it is conserved 
only in the mean. This model presents a first order tran- 
sition for all values of p except for p = 0.5, for which 
the phase transition is second order. In Ref. [Of we con- 
sidered a canonical framework for this model, the Ising 
model with fixed magnetization (IMFM). This constraint 
is important because it allows a direct exploration of the 
coexistence region, which is forbidden in the grand canon- 
ical formalism, but physically relevant because the sys- 
tem seen as an ensemble of interacting particles can take 
values of (T, p) in this region. Since we are interested in 
finite crossovers, the ensemble in which the theoretical 
model is formulated is important when we try to match 
it with experimental results. 

In the IMFM the number A of particles which are lo- 
cated on a lattice is strictly fixed in a fixed finite volume 
V — L^. The Hamiltonian reads 



H, 



IMFM 



Pt 



A 

i=l 






(1) 



±1}. 



where Vq is a constant potential strength and {ai 
The interaction acts between nearest neighbours. 

The Hamiltonian can be rewritten in terms of Si — 
{ai + f )/2, the total number of particles is ^^ Si = A, 
and the density p = A/V is fixed as a constraint so that 
the partition function in the canonical ensemble reads 
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where /3 = T~^ is the inverse temperature. 

In Ref. iQ we worked out this model in Sd for finite 
systems with hnear dimensions ranging from i = 10 to 
L = 48. The energy samphng of simulation events and 
the determination of critical exponents relying on finite 
size scaling (FSS) assumptions |lj,|l^ lead to the conclu- 
sion that the system experiences a second order transition 
at every value of the density p. 

This result was in contrast with the common under- 
standing that the LGM has a first order transition line 
except for a single point located sX p — 0.5. The main 
argument to explain this difference was the following. If 
one looks at the phase diagram and lowers the temper- 
ature for fixed density from above the transition line in 
the canonical description of this model it is possible to 
cross this line and enter into the 'coexistence' zone. In 
the grand canonical formulation (LGM) |Q, the states 
of this zone are not accessible as equilibrium states. In- 
deed, in the IMFM the magnetisation m — ^,- Ui which 
is equivalent to the density p is not discontinuous at the 
separation line as it should be in a first order transition. 
However, we could not firmly establish the nature of the 
transition in the thermodynamic limit because of limi- 
tations imposed by numerics on the size of the largest 
system which could be generated. This is the reason for 
which we discuss different tests aimed to characterise the 
smoothed crossover. 

The observed behaviour of the caloric curve and the 
specific heat cannot exclude a possible non-homogeneity 
in the 'coexistence' zone and the importance of surface 
effects. Hence our first test concerns the topology of the 
finite system in the vicinity of the line which separates 
the two phases. Fig. |l| shows two configurations gener- 
ated in the framework of the IMFM for a 2d system with 
linear size L = 400. They correspond to two values of the 
density, just below the transition line (low temperature 
side). In the ordinary Ising model, the case p = 0.5 cor- 
responds to a second order phase transition, and p = 0.3 
to a first order transition. In the second case, the non- 
homogeneous configurations arc those which dominate in 
the thermodynamic limit below the transition. One ex- 
pects the system to be divided in two well-defined phases 
separated by a transition domain of linear dimension ^, 
where ^ is the correlation length. However, as one can see 
in Fig. |l|, in the IMFM the system looks rather homoge- 
neous in both cases and one expects a self-similar pattern 
at different scales which is a qualitative indication of a 
continuous behaviour. At very low temperatures, the sys- 
tem is made up of large and compact clusters and hence 
gets non-homogeneous, but we checked that this happens 
for p ~ 0.5 as well as for any other value of p. 

The next point concerns the caloric curve for the 
IMFM in id. In the thermodynamic limit the caloric 
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FIG. 1. Distribution of particles in a 2d system calculated 
in the framework of the canonical IMFM, at the phase sepa- 
ration line for densities (a) p — 0.5 and (b) p = 0.3. The dark 
areas correspond to particles, the grey areas to vacuum. 



curve of a system which undergoes a first order tran- 
sition shows a plateau which signals the generation of 
latent heat. This is why the first experimental caloric 
curves [El, apparently showing a plateau, were considered 
as a reminiscence of a liquid-gas phase transition present 
in nuclear matter. In the framework of the IMFM, we 
constructed the caloric curve for different sizes of 3d sys- 
tems. Fig. g shows these curves for L = 10, 24 and 48 and 
a density p = 0.3. All three curves exhibit an inflection 
point and the slope of the curve in the interval of energy 
where it increases rapidly gets steeper with increasing L. 
However, one observes only a very small change between 
L = 24 and L = 48 which seems to indicate that the 
asymptotic limit is close. It is nevertheless clear that 
for i = 48 there is no sign for the appearance of a real 
plateau. Thus for this system one observes a behaviour 
of the caloric curve which indicates a continuous transi- 
tion between two homogeneous phases in the asymptotic 
limit. Such a behaviour is in qualitative agreement with 
recent data analyses Q. Similarly to the FSS analysis, 
it again does not allow to establish what happens in the 
infinite system. 

As a third test we introduce the microcanonical ap- 
proach in order to compare the caloric curve with the 
same curve obtained in the framework of the canonical 
ensemble. It is well established that a closed thermalised 
system with fixed energy and number of particles, which 
is a finite microcanonical ensemble shows a caloric curve 
which is multivalued in energy for fixed temperature, i.e. 
shows an oscillation ("S" curve behaviour). As a conse- 
quence the specific heat gets negative for fixed volume 
over certain intervals of energy ||l6|,|l^ • This effect is due 
to the non- homogeneity of the system which characterizes 
a first order transition. The surface energy of the clusters 
in the coexistence phase can be read from the area of the 
domains obtained by means of a Maxwell construction in 
the region of energy where the two phases coexist [Q. 
The canonical caloric curve does not show this effect. 
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FIG. 2. Canonical caloric curves for 3d IMFM systems of 
linear sizes L — 10, 24, 48 with periodic boundary conditions 
and a density p = 0.3. 



In order to test the existence or absence of this phe- 
nomenon in our model we introduce a Metropolis Monte 
Carlo microcanonical algorithm |]18[. We apply it to 
the g-state Potts model [|l9[ in two dimensions, which 
presents a second order transition for q < 4 and a first 
order transition for q > 5 and is used here as a test 
model. As expected, the system shows the two types of 
caloric curves, i.e. the characteristic "S" behaviour of the 
caloric curve when the transition is first order [q — 5), 
and a monotonous rise when the transition is second or- 
der ((7 = 4), see Fig. ^. This can be compared to a canon- 
ical simulation of the Potts model which is also shown in 
the figure. The backbending is absent for g = 5. 

In the case of the 3d IMFM, precise canonical and 
microcanonical simulations lead to the results shown in 
Fig. ^. The figures correspond to densities of p ~ 0.3 and 
p — 0.5. There exists no sign for the existence of multival- 
uedness of the energy for fixed temperature which would 
indicate a first order transition. We also observe from 
these figures that microcanonical and canonical results 
are already the same for a still relatively small lattice 
volume (L — 10). This result supports the previous indi- 
cations of a continous transition in the thermodynamic 
limit. 

These results raise the general question of the role 
played by constraints on the order of the transition. In 
order to investigate this point further we consider the 
Potts model for g > 5 and introduce different constraints. 
In practice we fix the number N{p) of spins in a given 
direction (p = 0, . . . , g — 1) to a given value and gener- 
ate the caloric curve for the systems in the framework of 
the microcanonical ensemble. Some results are shown in 
Fig. ||. One sees that the curves no longer correspond to 
a multivalued temperature T as a function of energy but 
show a smooth increase, with the presence of an inflec- 
tion point as expected for a second order transition when 
one goes to the thermodynamic limit. This result goes 
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FIG. 3. Caloric curves for the q — 4 and 5 = 5 Potts model 
in a 16 X 16 system. Black points (full line) are obtained from 
a microcanonical simulation, white points (dotted line) from 
a canonical simulation. See comments in the text. 



the same direction as the result obtained from the IMFM 
and confirms indeed the fact that constraints may have 
drastic effects on the order of the transition. 

In summary, different tests show that the IMFM, in 
its microcanonical or canonical formulation, exhibits the 
characteristic features of an homogeneous system which 
will undergo a smooth transition from one phase to the 
other in the thermodynamic limit. It may be surpris- 
ing to find such a behaviour in the framework of the 
LGM, but the constraints imposed on the model seem 
to produce this effect. The strongest argument point- 
ing towards this conclusion concerns the microcanonical 
test which indicates that the expected effect character- 
izing a discontinuous transition on the caloric curve is 
not seen. This is also in agreement with the finite size 
scaling analysis carried out in Ref. [|3|. It shows that 
microcanonical calculations may be of help in the char- 
acterisation of the phase crossover in finite systems. As 
already stated above, we however cannot conclude that 
the phase transition is indeed necessarily of second or- 
der in the thermodynamic limit. Hence it may also be 
difficult to decide about the order of a phase transition 
in infinite nuclear matter from experimental results ex- 
tracted from finite nuclei collisions. 

The sharpest experimental test would correspond to 
measurements in which energy and number of particles 
are strictly fixed and hence the fragmenting system could 
be assimilated to a microcanonical ensemble. Its caloric 
curve would then be able to reveal not only the charac- 
teristics of the crossover in the finite system but also the 
nature of the thermodynamic transition. Very recently 
it was claimed that the analysis of fragmentation events 
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FIG. 4. Caloric curve of the d = 3 IMFM with L = 10 at 
p = 0.3 (upper part) and p = 0.5 (lower part). Canonical 
and microcanonical results are undiscernable. See comments 
in the text. 

show negative values for the specific heat pO] . The confir- 
mation of these results would be a proof for the presence 
of a first order phase transition and confirm the predic- 
tions of models which where proposed in the near past 
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